ON  THE  GLOBAL  REGULARITY  OF 
SUB-CRITICAL  EULER-POISSON  EQUATIONS  WITH  PRESSURE 

EITAN  TADMOR  AND  DONGMING  WEI 


Abstract.  We  prove  that  the  one-dimensional  Euler-Poisson  system  driven  by  the  Poisson  forcing 
together  with  the  usual  7-law  pressure,  7  >  1,  admits  global  solutions  for  a  large  class  of  initial  data. 
Thus,  the  Poisson  forcing  regularizes  the  generic  finite-time  breakdown  in  the  2x2  p-system.  Global 
regularity  is  shown  to  depend  on  whether  or  not  the  initial  configuration  of  the  Riemann  invariants 
and  density  crosses  an  intrinsic  critical  threshold. 


1.  Introduction 


It  is  well  known  that  the  systems  of  Euler  equations  for  compressible  flows  can  and  will  breakdown 
at  a  finite  time  even  if  the  initial  data  are  smooth.  A  prototype  example  for  such  systems  is  provided 
by  the  2x2  system  of  isentropic  gas  dynamics 


(1.1) 


Pt  +  ( pu)x  =  0 
(pu)t  +  (pu2) x  =  -px  , 


where  the  pressure  p  =  p(p)  is  given  by  the  usual  7-law,  p(p)  =  Ap1 .  By  using  the  method  introduced 
in  [La64]  to  deal  with  pairs  of  conservation  laws,  it  can  be  shown  that  (1.1)  will  lose  the  C1  -smoothness 
due  to  the  appearance  of  shock  discontinuities  unless  its  two  Riemann  invariants  are  nondecreasing. 
Thus,  the  finite  time  breakdown  of  (1.1)  is  generic  in  the  sense  that  it  holds  for  all  but  a  “small  set” 
of  initial  data. 

On  the  other  hand,  if  we  replace  the  pressure  by  Poisson  forcing,  then  we  arrive  at  the  system  of 
Euler-Poisson  equations 


(1.2) 


Pt  "P  ( pu)x  —  0  , 

(pu)t  +  (pu2) x  =  —kppx  k>  0  , 


subject  to  initial  data  (uq  ,  po  >  0) .  Here  ip  =  <p(p)  is  the  potential,  which  is  dictated  by  the  (one¬ 
dimensional)  Poisson  equation,  <pxx  =  —p .  In  this  case,  there  is  a  “large  set”  of  initial  configurations 
which  yield  global  smooth  solutions.  More  precisely,  [ELT01]  have  shown  that  (1.2)  admits  a  global 
smooth  solution  if  and  only  if 


(1.3) 


Uqx(x)  >  -y/2kp0(x)  . 


Thus,  following  the  terminology  of  [LT02],  the  curve  uqx  +  \/2kf)Q  =  0  is  a  “critical  threshold”  in 
configuration  space  which  separates  between  initial  configurations  leading  to  finite  time  breakdown 
and  a  “large  set”  of  sub-critical  initial  configurations  which  yield  global  smooth  solutions.  In  particular, 
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(1.3)  allows  negative  velocity  gradients  (-  depending  on  the  local  amplitude  of  the  density),  which 
otherwise  are  excluded  in  the  case  of  inviscid  Burgers  equations,  corresponding  to  k  =  0. 

In  this  paper  we  turn  our  attention  to  the  full  Euler-Poisson  equations  driven  by  both  -  pressure 
and  Poisson  forcing, 


(1.4) 


Pt  +  {pu) x  =  0  , 

<  (pu)t  +  {pu2) x  =  -p{p) x  -  kppx,  k>0, 
PxX  —  P  ■ 


These  equations  govern  different  phenomena,  ranging  from  the  largest  scale  of  e.g.,  the  evolution 
gravitational  collapse  in  stars,  to  applications  in  the  smallest  scale  of  e.g.,  semi-conductors.  There  is 
a  considerable  amount  of  literature  available  on  the  local  and  global  behavior  of  Euler-Poisson  and 
related  problems.  Consult  [Ma86]  for  local  existence  in  the  small  7Ls-neighborhood  of  a  steady  state 
of  self-gravitating  stars,  [CW96]  for  global  existence  of  weak  solutions  with  geometrical  symmetry, 
[Gu98]  for  global  existence  for  3-D  irrotational  flow,  [MN95]  for  isentropic  case,  and  [JR00]  [PRV95] 
for  isothermal  case.  Consult  [Pe90]  [MP90]  [Si85]  [WC98],  [BW98]  for  non-existence  results  and 
singularity  formation.  The  question  of  global  smoothness  vs.  finite  breakdown  was  studied  in  recent 
series  of  works  of  Liu  and  Tadrnor,  in  terms  of  a  critical  threshold  phenomena  for  1-D  “pressure-less” 
Euler-Poisson  equations,  [ELT01]  and  2-D  restricted  Euler-Poisson  equations,  [LT02,  LT03]. 

The  natural  question  that  arises  in  the  present  context  of  full  Euler-Poisson  equations  (1.4)  is 
whether  the  pressure  enforces  a  generic  finite  time  breakdown  or,  whether  the  presence  of  Poisson 
forcing  preserves  global  regularity  for  a  “large”  set  of  initial  configurations.  We  answer  this  question 
of  “competition”  between  pressure  and  Poisson  forcing,  proving  that  the  Euler-Poisson  equations  (1.4) 
with  7  >  1  admit  global  smooth  solutions  for  a  “large  set”  of  sub-critical  initial  data  such  that 

(1-5)  u0x(x)  >  -ilovVo(x)  +  ,  7>1. 

p0{x)  2 

Here,  Kq  is  a  constant  depending  on  k,  7  and  the  initial  data.  In  the  particular  (and  important)  case 
of  isothermal  equations,  7  =  1,  we  have  Kq  =  \/2k  and  (1.5)  amounts  to  a  sharp  critical  threshold, 

(1.6)  u0x(x)  >  ~^2kp0(x)  +  y/A^POxj-X^ ,  7=1. 

Po{x) 

The  inequalities  (1.5), (1.6)  quantify  the  competition  between  the  destabilizing  pressure  effects,  as  the 
range  of  sub-critical  initial  configurations  shrinks  with  the  growth  of  the  amplitude  of  the  pressure, 
A,  while  the  stabilizing  effect  of  the  Poisson  forcing  increases  the  sub-critical  range  with  a  growing  k. 
In  particular,  (1.6)  with  A  =  0  recovers  the  pressure-free  critical  threshold  (1.3). 

The  paper  is  organized  as  follows.  In  section  2,  we  reformulate  the  system  (1.4)  with  its  Riemann 
invariants  as  a  preparation  for  the  analysis  carried  out  in  sections  3  and  4.  In  section  3,  we  prove 
our  main  results,  providing  sufficient  conditions  for  “large  sets”  of  sub-critical  initial  configurations 
which  yield  global  smooth  solution.  In  section  4,  we  give  examples  of  finite  time  breakdown  for  super¬ 
critical  initial  data.  Combining  our  results  in  sections  3  and  4,  they  confirm  the  existence  of  a  critical 
threshold  phenomena  for  the  full  Euler-Poisson  equations  (1.4). 
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2.  Riemann  Invariants 


2.1.  The  Euler-Poisson  equations  with  7-law  pressure:  7  >  1.  We  begin  by  rewriting  the 
Euler-Poisson  equations  (1.4)  as  a  first  order  quasilinear  system 


where  the  Jacobian  J  :  = 


A'yp 7  2  u 


has  two  different  eigenvalues 


A  :=  u  —  \J~Arfp:L^~  <  p  :=  u  +  y 1  A'yp2^-  . 

and  let  R  and  S  denote  the  Riemann  invariants  of  the  corresponding  Euler  system  (1.1) 

.  .  2v/ A'y  1 zl  ,  „  2  + A'y  t-i 

(2.2)  R:=u - p  2  and  S:=u-\ - p  2  . 

7  —  1  7  —  1 

They  satisfy  the  coupled  system  of  equations, 

(2.3a)  Rf  +  A Rx  =  —k(px  , 

(2.3b)  St  +  nSx  =  -kipx , 

coupled  through  the  Poisson  equation  —(j)xx  =  p .  If  we  set  r  :=  Rx  ,  s  :=  Sx  then  upon  differentiation 
of  (2.3)  we  get 


(2.4a) 

(2.4b) 


n  +  A rx  +  A srs  +  A Rr2  =  kp  , 
st  +  l^sx  +  pss 2  +  pRrs  =  kp  . 


Next,  we  observe  that  A  =  -^5  _  —  R)  and  p  =  ^±±1  7.  3_1  (5  _  Ry  Hence,  expressed  in  terms 

7  —  1 

of  6  :=  — — — ,  we  have  for  7  >  1, 

X  1  +  0  j  \  1-0  7  —  1  ^  n 

*R  =  HS  =  —y—  and  XS  =  pR  =  — - — ,  6  :=  — —  >  0, 

and  the  pair  of  equations  (2.4)  is  recast  into  the  form 


(2.5a) 

(2.5b) 


v  1  +  9  2  1  —  d 

r  H - - — r  H - - — rs  =  kp , 

,  1  +  0  2  1-0 
s  H - - — s  H - - — rs  =  kp  . 


Here  and  below  {}'  :=  dt  +  \dx  and  j}7  :=  dt  +  pdx  denote  differentiation  along  the  A  and  p,  particle 
paths, 

Ta  :=  {(x,t)  |  x(t)  =  \(p(x,t),u(x,t))},  :=  {(x,t)  \  x(t)  =  p(p(x,t),u(x,t))}. 

To  continue,  we  rewrite  the  equation  for  p  as 

(2-6)  (pt  +  A px)  +  — - — px  +  p — — —  =  0  . 

( , .  /— —  ^—3  ,  p  —  A  r— —  ^  — 1  S  —  V 

bmce  s  —  r  =  dx  —  Rx  =  2^  A'yp  2  px  ,  it  enables  us  to  express  — - — px  =  y  A'yp  2  px  =  p — - —  ,  so 

that  the  p  equation  (2.6)  can  be  written  along  the  A  particle  path  as  ps  +  ps  =  0.  Similarly,  it  can  be 


4 


EITAN  TADMOR  AND  DONGMING  WEI 


written  along  the  //  particle  path  as  p'  +  pr  =  0.  Assembling  the  above  equations  together,  we  arrive 
at  the  following  system  governing  r,  s  and  p , 


(2.7a) 


\  1  + 

r  H - 


■r  + 


1  - 


-rs 


2 

p'  +  ps 


=  kp, 
=  0, 


and 


(2.7b) 


/  1  + 
s  H - 


+ 


1  - 


-?’S  = 


y  p'  +  pr  = 

Finally,  we  use  the  integration  factors  l/y/p  and  r/2pv/p  in  the 
in  (2.7),  to  conclude 


kp , 

0. 

first  and  second  equations  of  each  pair 


(2.8a) 

(2.8b) 


(  r  V  1  +  6  r2  9  rs 
(  s  \'  1  +  9  s2  9  rs 

\^p)  ^  2  ~  T> ' 


2.2.  The  isothermal  case  7  =  1.  In  this  case,  the  two  eigenvalues  are  A  =  u  —  \fA  <  p  =  u  +  \J~A 
with  the  corresponding  Riemann  invariants  R  =  u  —  \[A In  p  and  S  =  u  +  \Z~A In  p.  Their  derivatives, 
r  and  s,  satisfy  the  pair  of  equations,  corresponding  to  (2. 8a), (2. 8b)  with  9  =  (7—  l)/2  =  0, 


(2.9a) 


(2.9b) 


1  r 

Is2 

+  2~P  =  k^- 


3.  Global  smooth  solutions  for  sub  critical  initial  data 

For  the  pressure-less  Euler-Poisson  equations  (1.2),  the  evolution  of  ux  and  p  could  be  traced 
backwards  along  the  same  particle  path  to  their  initial  data  at  t  =  0.  The  scenario  becomes  more 
complicated  with  the  additional  pressure  term,  due  to  the  coupling  of  r  and  s  along  different  particle 
paths  which  are  traced  back  to  different  neighborhoods  of  the  initial  line  t  =  0.  This  is  the  main 
obstacle  in  finding  the  sharp  critical  threshold  of  the  full  Euler-Poisson  system  (1.4).  To  this  end,  we 
will  seek  invariant  regions  for  the  coupled  system,  governing  the  Riemann  invariants.  We  begin  this 
section  with  the  following  lemma. 

Lemma  3.1.  Given  that  the  total  charge  Eq  :=  po(x)dx  <  00,  then  p(x,t)  and  u(x,t)  remain 
uniformly  bounded  for  all  t  >  0. 

Proof.  Under  the  given  condition,  we  can  set  (e.g.,  [ELT01,  p.  116]) 

tPx(x,t)  =  ^  J  p(£,  t)d£,  —  J  p(£,  t)d^j  , 

which  satisfies  —E0  <  px(x,  t)  <  Eq  ,  for  all  f  >  0  and  x  €  M . 

Recall  the  transport  equations  (2. 3a), (2. 3b)  which  govern  the  Riemann  invariants  along  different  char¬ 
acteristics  R'-\-kpx  =  S'  +  kipx  =  0.  Since  <px  is  bounded,  these  transport  equations  tell  us  that  R  and 
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S  remain  uniformly  bounded  with  at  most  a  linear  growth  in  time.  Indeed,  for  all  M  3>  1  we  have 
(3.1)  sup  {\R(x,t)\,\S(x,t)\}<C0  +  kE0t,  C0  :=  sup  {|-R0(aOI,  |50(x)|} . 

\x\ <M  \x\<M -\-Uoot 

Take  the  sum  and  difference  of  S  and  R  to  find  that  u(x,  t )  and  p(x,  t )  in  (2.2)  remain  bounded, 


(3.2)  tioo  :=  sup  |tt(x,  i)|  <  Co  +  kEot, 

\x\<M 


sup  p(x,  t )  <  Const. 

\x\<M 


( Co  +  kEot )~,  7  >  1, 

exp  (kEot).  7=1. 

□ 


We  note  in  passing  that  the  time  growth  asserted  in  (3.2)  is  probably  not  sharp;  the  estimate  can  be 
improved  after  taking  into  account  the  uniform  bounds  of  Rx/yfp  and  Sx/ yfp  discussed  in  theorems 
3.1  and  3.2  below. 


Remark  3.1.  According  to  Lemma  3.1,  the  only  way  that  the  full  Euler-Poisson  system  (1-4)  breaks 
down  at  a  finite  time  is  through  the  formation  of  shock  discontinuities  where  \ux\  and/or  \px\  blow  up 
|  oo ,  but  neither  will  concentrate  at  any  critical  point.  This  is  in  contrast  to  the  breakdown  of  the 
“pressure-less”  Euler-Poisson  equations  (1.2),  where  —ux(x,t)  =  p(x,t)  |  oo  simultaneously  at  the 
critical  time. 


3.1.  Critical  Threshold  for  isothermal  case:  7  =  1.  We  begin  with  the  isothermal  case,  7  =  1, 
which  plays  an  important  role  in  various  applications.  Compared  with  the  general  case  (2.8),  the 
isothermal  case  becomes  simpler  due  to  the  fact  that  9  =  0  decouples  the  dependence  on  r  and 
s  through  the  mixed  term  9rs,  which  disappears  from  left  hand  side  of  (2.9).  Here  we  prove  the 
following  sharp  characterization  of  the  critical  threshold  phenomena. 

Theorem  3.1.  Consider  the  isothermal  Euler-Poisson  system  (1-4)  with  pressure  forcing  p(p)  =  Ap , 
and  subject  to  initial  data  (uq,  p$  >  0)  with  finite  total  charge,  Eq  =  po(x)dx  <  00.  The  system 
admits  a  global  smooth,  C1  -solution  if  and  only  if 

(3.3)  u0x(x )  >  —  \J 2kpo(x)  +  ,  ViGl. 

Po(x) 

Remark  3.2.  Expressed  in  terms  of  the  Riemann  invariants  specified  in  §2.2,  ux  ±  \/Apx/p,  theorem 
3.1  states  that  the  isothermal  Euler-Poisson  equations  admit  global  smooth  solutions  for  sub-critical 
initial  conditions, 

(3.4)  so  >  —\jTkp§  and  ?’o  >  —\r2kpo. 

T  S 

Proof.  We  define  X  := -  and  Y  :=  — .  Equations  (2. 9a), (2. 9b)  then  read 

s/p  sfp 


(3.5a) 

X'  =  i^-iflk-X2), 

(3.5b) 

<N 

1 

CM 

Tf- 

II 

It  follows  that 

>0,  Xe(-s/2k,V2k) 

X'  < 

=  0,  \X\  =  y/2 k, 

<0,  |  X\>/2k, 
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and  similarly, 

"  >  0,  Ye  (~V2k,V2k), 

Y'  l  =0,  \Y\  =  y/Zk, 

<  0  ,  |Y|  >  V2k  . 

Thus,  starting  with  (3.4),  W0,  Yq  >  —  \[2 ~k,  we  find  that  X  and  Y  remain  bounded  within  the  invariant 
region  [~y/2k,y/2k],  or  otherwise,  they  are  decreasing  outside  this  interval.  We  conclude  that 

X(-,t),Y(-,t)  <  maxjv^XoO),  Y0(-)}  • 

Lemma  3.1  tells  us  that  p  is  bounded.  The  boundedness  of  X ,  Y  and  p  imply  that  r  =  X^fp  and 
s  =  Y \J~p  remain  bounded  for  all  t  <  oo,  and  hence  the  Euler-Poisson  system  (1.4)  admits  a  global 
smooth  (^-solution. 

Conversely,  suppose  that  there  exists  X0  =  X(xq)  <  —\[2 ~k.  We  will  show  that  this  value  will  evolve 
along  rA(x0,  0)  such  that  X(-,t)  will  tend  to  — oo  at  a  finite  time.  To  this  end,  assume  that  Y  is  well 
behaved,  i.e. ,  Lq(-)  >  —  \[2 \k  so  that  Y(-,t)  <  Y\  :=  max  jlo(-),  \/2A;|  for  all  Vs  (otherwise,  the  finite 
time  blow  up  of  Y  can  be  argued  along  the  same  lines).  It  follows  that  s  =  Y s/p  <  Y\^fp  and  inserting 
this  into  p'  =  —  ps,  we  find  ps  >  —  Yip3^2.  This  yields  the  lower-bound 

/  2  N  2 

P  > 


Y\t  +  2/yfpQ, 

and  together  with  (3.5a),  we  conclude  that  X(-,t)  satisfies  the  following  Ricatti  equation  along  the 
TA-path, 

X\ 


Y\t  +  2/ ^/po 


Xz 


X\  :=  (X0J  -  2k) /X^  >  0. 


(3.6)  Xv  <  - 

Integration  of  (3.6)  yields 

“ V  ’  V  "  Xx  In  (1  +  v^Yit/2)  +  YxX 0 

Thus,  starting  with  Xq  <  —\/2k  <  0  it  follows  that  there  exists  a  finite  critical  time  tc  >  0  such  that 
X(t  |  tc)  tends  to  — oo.  □ 


Y 


The  critical  threshold  condition  (3.3)  reflects  the  competition  between  the  Poisson  forcing  and  the 
pressure.  It  yields  global  smooth  solutions  for  a  “large”  set  of  initial  configurations  allowing  negative 
velocity  gradients.  In  the  particular  case  that  there  is  no  pressure,  A  =  0  ,  (3.3)  is  reduced  to  the  critical 
threshold  condition  of  the  “pressure-less”  Euler-Poisson  equations  uqx  >  —y/,2kpo(x)  of  [ELT01]. 

3.2.  Critical  threshold  for  7  >  1.  The  equations  for  the  Riemann  invariants  (2.8a),  (2.8b)  are 
coupled  through  the  mixed  term,  Qrs/2.  We  note  in  passing  that  it  is  possible  to  get  rid  of  this  mixed 
term  when  integrating  (2.7a),  (2.7b)  with  the  integration  factors  p^7-3^4,  and  ?’p^7_'^4( 3  —  y)/4  in 
the  first  and  second  equations  in  each  pair,  yielding 


1  0  2  1 

+  — r  P  2 

1  T  0  2  1 

+  —  s» 1 


7 

kp  2  . 
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Nevertheless,  it  will  prove  useful  to  use  the  same  integration  factors,  1  / ^/~p  and  r^p^fip  which  led 
to  (2.8).  The  main  task  is  to  identify  the  invariant  region  associated  with  (2.8),  corresponding  to  the 
isothermal  invariant  region  [— v^2 ~k,  V2k]  discussed  in  theorem  3.1. 

Theorem  3.2.  Consider  the  Euler-Poisson  system  (1-4)  with  7  law  pressure  p(p)  =  Ap1 ,  7  >  1, 
subject  to  initial  data  (uq,  po  >  0)  with  finite  total  charge,  Eq  =  po(x)dx  <  00.  Then,  there  exists 
a  constant  Kq  >  0  depending  on  k,  7  and  the  initial  conditions  (specified  in  (3.9b)  below),  such  that 
the  Euler-Poisson  equations  (1-4)  admit  a  global  smooth,  C1  -solution  if, 

\Pox(x)\ 


(3.8) 


u0x(x)  >  -K0y/p0(x)  +  y/Ay 


Po(x)32J 


Before  we  turn  to  the  proof  of  this  theorem,  several  remarks  are  in  order. 

Remark  3.3.  Expressed  in  terms  of  the  Riemann  invariants,  r  =  ux  —  y/~AypQXf Pq  7^2  and  s  = 
ux  +  '/A^Pox/p<o  7'*//2;  the  critical  threshold  (3.8)  reads 
n  \  ro(x )  s0(x) 

[  j  VMx)'  VMx) 

The  constant  Kq  is  given  by 

—6  Mq  +  y/WMf+mTTe) 


>  —Kq 


(3.9b) 


Krt  = 


2(1  +  0) 


M0  =  max  <  V2k. 


r  o(x)  s0(x) 


VPo(x)  vw) 

We  mention  two  simplifications  which  are  summarized  in  the  following  two  corollaries.  We  first 
observe  that  if  the  initial  configurations  satisfy  the  upper-bound  r0(x),  sq(x)  <  y/2kpo(x)  then  (3.9b) 

yields  Mq  =  \[2h,  hence  Ko  =  ——7;,  and  theorem  3.2  implies  the  following. 

1  +  0 

Corollary  3.1.  Consider  the  Euler-Poisson  system  (1-4)  with  7  law  pressure  p(p)  =  Ap1 ,  7  >  1, 
subject  to  initial  data  (uq,  po  >  0)  with  finite  total  charge,  Eq  =  po{x)dx  <  00.  Then,  the  Euler- 
Poisson  equations  (1-4)  admit  a  global  smooth,  C1  -solution  if  for  all  x  G  R, 


(3.10) 


\uox(x) I  <  y/2kp0(x)  -  y/Ay 

Po{x 


-7 

2 


0MO  -  (0MO  +  y/8k(l  +  0))/\/2 


The  next  result  follows  from  the  trivial  inequality  —Kq  <  , 

H  J  u“  2(1  +  0) 

Corollary  3.2.  Consider  the  Euler-Poisson  system  (1-4)  with  a  y-law  pressure  p(p)  =  Ap 7  ,  7  >  1, 
subject  to  initial  data  (uq,pq  >  0)  with  finite  total  charge,  Eq  =  J//  po(x)dx  <  00.  Then,  the  Euler- 
Poisson  equations  (1-4)  admit  a  global  smooth,  C1  -solution,  if  for  all  x  G  R, 


u0x(x)  >  - 


(3.11) 


7-1 

,  7 - -  max 

x/2^  2(7+1)  X 


y/2kpo(x),  U0x(x)  +  y/Ay 

Po(x 


-7 

2 


\P0x(x)\ 
Po{x ) V 


Remark  3.4.  We  observe  that  as  in  the  isothermal  case,  the  critical  threshold  in  its  various  versions 
(3. 8), (3. 9),  (3.10)  and  (3.11),  allow  a  “large  set”  of  initial  configurations  with  negative  velocity  gradi¬ 
ent,  due  to  the  competition  between  the  stabilizing  Poisson  forcing  kpfi(p)x  and  the  destabilizing  pres¬ 
sure  A(p'r)x.  In  the  extreme  case  that  Poisson  forcing  is  missing  k  =  0 ,  the  breakdown  of  the  system 
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is  generic  unless  uqx  is  positive  enough  (so  that  cq,  so  >  Oj.  In  the  other  extreme  of  a  “pressure-less” 
Euler- Poisson,  A  =  0, 7=1,  the  critical  thresholds  (3.8),  (3.10)  are  reduced  to  uqx(x)  >  —  yi/2kpo(x) , 
which  coincides  with  the  “pressure-less”  critical  threshold  (1.3)  found  in  [ELT01]. 

T  S 

Proof.  Expressed  in  terms  of  X  :=  -  and  Y  :=  — — ,  equations  (2.8)  read 

Vp  Vp 

(3.12a)  X'  =  J-p(k  -  ^X2  +  , 

(3.12b)  y'  =  jp{k  -  1-^-y2  +  d-XY^j  . 

We  seek  an  invariant  region  of  the  form  [—Kq,  Mq\,  with  Kq,  Mq  >  0  yet  to  be  determined.  We  begin 
by  noticing  that  if  X,Y  <  M  then1  X+Y  <  M2,  and  recalling  that  6  >  0,  (3.12)  then  yield 

<  y/p(k  ~  l~^X2  +  -2M2)  ■  X>0, 

y  <  Vp(*  -  X-^-Y2  +  e-M2^J  ,  Y  >  0. 

This  in  turn  implies  that 

(3.13)  X  and  Y  are  decreasing  if  X,  Y  >  C+,  C+  =  C+(M)  := 

The  solution  of  C+{M )  =  M  yields  M  =  \Pl \k.  Thus,  X  and  Y  are  decreasing  whenever  X,  Y  >  M  =  \/~Yk, 
and  we  end  up  with  the  upper-bound 


(3.14) 


X(-,t),Y(-,t)  <  Mq,  Mq  :=max|\/2 k,  Xq(x),Yq(x)^  . 


In  a  similar  manner,  we  study  the  lower  bound  of  the  invariant  region.  By  (3.14)  and  (3.12)  yield 


(3.15a) 


(3.15b) 


X'  >  Vp(k  ~  l-^~X2  +  \mqX\  ,  X  <  0, 
y  >  V p(k  -  ^ Y 2  +  6-MqY^  ,  Y  <  0, 


which  in  turn,  imply  that 

(3.16a)  X  and  Y  are  increasing  if  0  >  X,  Y  >  —Kq, 

where  Kq  is  the  smallest  root  of  the  quadratics  on  the  right  of  (3.15), 


(3.16b) 


Kq  ■.= 


2(1  +  0) 


The  critical  threshold  condition  (3.8)  tells  us  that  at  t  =  0,  Xq,Yq  >  —Kq  and  (3.16a)  implies  that 
X(-,t)  and  Y(-,t)  remain  above  the  same  lower-bound,  (3.8).  As  before,  the  bounds  of  X,  Y  and  p 
imply  that  r  =  X^fp  and  s  =  Y^fp  remain  bounded,  and  hence  the  Euler-Poisson  system  (1.4)  a 
global  smooth,  C'1-solution.  □ 


1We  let  Z+  =  max{X,  0}  and  Z-  =  min{Z,  0}  denote  the  positive  and  negative  part  of  Z . 
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4.  Finite  time  breakdown  for  super-critical  initial  data 


Consider  the  Euler-Poisson  system  (1.4)  with  a  7-law  pressure,  7  >  1,  and  subject  to  initial  data 
such  that  ?’o(x),so(x)  <  y/2k.  Then,  according  to  corollary  3.1,  the  following  critical  threshold  is 
sufficient  for  the  existence  of  global  smooth  solutions, 

uQx(x)  >  -\J 2kp0(x )  +  y/Ay-^^^-. 

Po{x)  2 

In  this  section  we  show  that  this  critical  threshold  is  also  necessary  for  global  regularity. 


Theorem  4.1.  Consider  the  Euler-Poisson  system  (1-4)  with  a  7 -law  pressure  p(p)  =  Ap1 ,  7  >  1, 
subject  to  initial  data  (tig,  po  >  0)  .  The  system  loses  the  C1  -smoothness  if  there  exists  an  x  6  M  such 
that 


(4.1) 


Uqx(x)  < 


-\f2kpo(x)  +  \fAy- 


\pQx(x)\ 
Po{x ) V 


Remark  4.1.  Expressed  in  terms  of  the  Riemann  invariants,  r  =  ux  —  \J Ay pox/ pg3  7^2  and  s  = 
ux  +  \J Ay pox/ Po  7^2,  the  condition  (4-1)  reads 

(4.2)  3x  G  M  s.t.  i’o(x)  <  —  y/2kpo (x) ,  or  so(x)  <  - y/2kpo (x)  . 

The  lack  of  smoothness  in  this  case  was  shown  in  theorem  3.1  for  7=1  and  is  extended  for  7  >  1 
below. 

T  S 

Proof.  Recall  equations  (3.12)  for  X  := -  and  Y  :=  — 

Vp  Vp 

(4.3a)  X'  =  J-p(k  -  ^ X 2  +  ^Xhj  , 

(4.3b)  Y'  =  Jp{k  -  +  ^Xbj  . 

In  the  proof  of  theorem  3.2,  we  have  shown  that  X  and  Y  have  an  upper  bound 


V2k,X0{x),Y0{x)y 

Suppose  that  there  exists  X0  =  X(xo)  <  —\/2k.  We  will  show  that  this  value  will  evolve  along 
Fa(xo,0)  such  that  X(-,f)  will  tend  to  —00  at  a  finite  time.  To  this  end,  assume  that  Y  is  well 
behaved,  i.e.,  YqO)  >  —\jTk  so  that  Y(-,t)  <  Mq  for  all  P s  (otherwise,  the  finite  time  blow  up  of  Y 
can  be  argued  along  the  same  lines).  It  follows  that  along  Ta^o,  0) 


(4.4) 


X(;t),Y(;t)<M0, 


Mo  :=  max 

X 


(4.5) 


X'  =  Jp(k  -  ^X2  +  9-XY^  <  y/p(k  -  -2X2^j  . 


Following  exactly  what  we  have  done  in  the  proof  of  theorem  3.1,  we  obtain  the  inequality 


(4.6) 


_ Mo _ 

Xi  In  (l  +  ^fpoMot/2)  +  MqXo  ’ 


where  X\  ■=  (Xg  —  2k) /Xq  >  0.  Thus,  starting  with  Xo  <  —\[2h  <  0  it  follows  that  there  exists  a 
finite  critical  time  tc  >  0  such  that  X(t  |  tc)  tends  to  —00.  □ 
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We  conclude  with  an  example  for  a  finite  time  breakdown. 
Example:  Suppose  at  t  =  0  ,  uq(x)  =  0  and 


Thus 


p0(x) 


1 ,  x  <  0  , 

1  -  Ye  ’  0  ^  x  ^  e> 


s0(x) 


-v^4 7(l-^)/2e,  0  <  x  <  e , 
0 ,  elsewhere . 


If  we  choose  e  small  enough,  then  so(x)  <  —  y/,2kpo(x)  for  0  <  x  <  e.  According  to  theorem  4.1,  the 
system  (1.4)  will  break  down  in  a  finite  time.  This  example  shows  that  even  if  the  fluid  is  near  rest 
at  t  =  0,  the  pressure  itself  could  still  lead  to  collision. 
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